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We analyze the stability of the Einstein static universe by considering homogeneous scalar pertur- 
bations in the context of f(R) modified theories of gravity. By considering specific forms of f(R), the 
stability regions of the solutions are parameterized by a linear equation of state parameter w = p/p. 
Contrary to classical general relativity, it is found that in f(R) gravity a stable Einstein cosmos with 
a positive cosmological constant does indeed exist. Thus, we are lead to conclude that, in principle, 
modifications in f(R) gravity stabilizes solutions which are unstable in general relativity. 
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I. INTRODUCTION 

Independent observations have confirmed that the Uni- 
verse is presently undergoing a phase of accelerated ex- 
pansion Although the introduction of a cosmolog- 
ical constant into the field equations seems to be the 
simplest theoretical approach to generate a phase of ac- 
celerated expansion, several alternative candidates have 
been proposed in the literature, ranging from dark energy 
models to modified theories of gravity. Amongst the lat- 
ter, models generalizing the Einstein- Hilbert action have 
been proposed. A nonlinear function of the curvature 
scalar, f(R), is introduced in the action given by 



S=— * / d 4 xV=g f(R) + S„ 



(1) 



where S m is the matter action. We consider k 2 = 8wG = 
1 throughout this work, for notational simplicity. Vary- 
ing the action with respect to g^ u provides the following 
field equation 

FR»u ~ \f9»u - V„ V„F + g^DF = T™ , (2) 

where F = df /dR. Note that the Ricci scalar is now a 
fully dynamical degree of freedom, which is transparent 
from the following relationship 



FR - 2f + 3 DF = T , 



(3) 



obtained from the contraction of the modified field equa- 
tion ([2]). One may generalize the action ([I]) by consid- 
ering an explicit coupling between an arbitrary function 
of the scalar curvature, R, and the Lagrangian density 
of matter 0]. Note that these couplings imply the vi- 
olation of the equivalence principle [3j, which is highly 
constrained by solar system tests. One may also mention 
alternative approaches, namely, the Palatini formalism 



'Electronic address: christian. bochmer@port.ac 
^Electronic address: lukas. hollenstein@port.ac. 
t Electronic address: francisco.lobo@port.ac.uk 



c.uk| 
ukj 



[U, [E| , where the metric and the connections are treated 
as separate variables; and the metric-affine formalism, 
where the matter part of the action now depends on the 
connection and is varied with respect to it |5j. 

A fundamental issue extensively addressed in the liter- 
ature is the viability of the proposed f(R) models [1,0]- 
In this context, it has been argued that most f(R) mod- 
els proposed so far in the metric formalism violate weak 
field solar system constraints [|[ , although viable models 
do exist @,M 03 [HI (Static and spherically symmetric 
solutions have also been found [Hj]). The issue of stability 
[HI also plays an imp ortant role in the viability of cosmo- 
logical solutions [TJ [HI, EH . In Ref. [h} it was argued 
that the sign of Jrr = d 2 f /dR 2 determines whether the 
theory approaches the general relativistic limit at high 
curvatures, and it was shown that for Jrh > the mod- 
els are, in fact, stable. The stability of the de Sitter so- 
lution in f(R) gravity has also been extensively analyzed 
in the literature [r|. Recently, f(R) models which have 
a viable cosmology were analyzed, and it was found that 
the models satisfying cosmological and local gravity con- 
straints are practically indistinguishable from the ACDM 
model, at least at the background level Note that 

to be a viable theory, the proposed model, in addition, to 
simultaneously account for the four distinct cosmological 
phases, namely, inflation, the radiation-dominated and 
matter-dominated epochs, and the late-time accelerated 
expansion [171 1 1 81 ] , should be consistent with cosmolog- 
ical structure formation observations pjjj. In the latter 
context, it has been argued that the inclusion of inhomo- 
geneities is necessary to distinguish between dark energy 
models and modified theories of gravity, and therefore, 
the evolution of density perturbations and the study of 
perturbation theory in f(R) gravity is of considerable 
importance [H [U, [Hj]. See Ref. [23] for studies of 
a parameterized growth factor approach to distinguish 
between modified gravity and dark energy models using 
weak lensing. 

In this work, we explore the stability of the Ein- 
stein static universe in f(R) modified theories of grav- 
ity. This is motivated by the possibility that the uni- 
verse might have started out in an asymptotically Ein- 
stein static state, in the inflationary universe context 
[24j . On the other hand, the Einstein cosmos has al- 
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ways been of great interest in various gravitational the- 
ories. In general relativity for instance, generalizations 
with non-constant pressure have been analyzed in [25[. 
In the context of brane world models the Einstein static 
universe was investigated in [261 ] . while its generalization 
within Einstein-Cartan theory can be found in (22|. Fi- 
nally, in the context of loop quantum cosmology, we refer 
the reader to [2§,|29j. 

We analyze the stability of the Einstein static universe 
against homogeneous scalar perturbations in the context 
of f(R) gravity. In the following section, we provide two 
specific forms of f(R), and analyze the stability of the 
solutions, by considering homogeneous scalar perturba- 
tions around the Einstein static universe. The stability 
regions are given in terms of the linear equation of state 
parameter w = p/p and the unperturbed energy density 
Po- 



ll. THE EINSTEIN STATIC UNIVERSE IN f(R) 
GRAVITY 



Consider the metric given by 
" dr 2 



-dr + a 2 {t) 



1 



r 2 {d6 2 + sin 2 



(4) 

For the Einstein static universe, a = clq = const, the 
Ricci scalar reduces to R = 6/a 2 ,, and the field equations 
take the following form 



Po 



f 
2 



2F 



Po 



f 
2 



(5) 



where po and po are the unperturbed energy density and 
isotropic pressure, respectively. 

In what follows, we consider specific forms of f(R), and 
analyze the stability against linear homogeneous scalar 
perturbations around the Einstein static universe given 
in Eqs. ©. Thus, we introduce perturbations in the 
energy density and the metric scale factor which only 
depend on time: 



p{t) = p (1 + Sp(t)) 



o(t) = a (1 + Sa(t)) . (6) 



Subsequently, we consider a linear equation of state, 
p(t) = wp(t), linearize the perturbed field equations and 
analyze the dynamics of the solutions. 



Firstly, motivated by the possibility that the universe 
might have started out in an asymptotically Einstein 
static state [HI], we analyze the case of f{R) oc R + R 2 , 
as in principle R 2 dominates for high curvatures, i.e. the 
early universe. Secondly, we consider the case of f(R) oc 
R + l/R, which is known to generate a late-time accel- 
erated expansion phase :3Q], and has been used in the 
weak-field limit constraints, as the l/R term dominates 
for low curvatures. 



A. f(R) oc R + R 2 theory 
Consider the case of 



f(R)=R- 



2„4 



act a 



o D 2 



(> 



R z - 2A . 



(7) 



where a = ±1, and a is a positive parameter. We intro- 
duced the factor <Zq/6 in the second term to considerably 
simplify the equations and results in the analysis outlined 
below. 

In this model the unperturbed field equations ([5]) take 
the form 



Po = -o + 3cra - A , 
% 

1 2 A 

PO = 2 + a0r + A ' 



(8) 
(9) 



and yield the following cosmological constant of the Ein- 
stein static universe 



A = ip (l + 3w) - 3aa 2 . 



(10) 



We derive an evolution equation for the scale fac- 
tor perturbation in the following way. The perturba- 
tions defined in Eq. ^ are introduced in the metric and 
the energy-momentum tensor. Then the perturbed field 
equations @ are linearized and the unperturbed equa- 
tions ijHJ) and © are subtracted to end up having only 
first order terms. The (tt)— component reduces to 

5p{t) = -3{l + w)5a(t) , (11) 

which we use to simplify the spatial component to 



[8aa 2 - p (l +w)(l + 3iu)] [-Acra 2 + p Q (l + w)] 2 Sa(t) 

- 2 [8aa 2 -p (l + w)] [Aaa 2 -p (l + w)] 6a" (t) + %oa 2 5a {A) (t) = . (12) 
I 

Note that in the general relativistic limit, a — > 0, Eq. (Tl2^) reduces to 

25a" (t) - p (l + w)(l + 3w)Sa(t) = , (13) 



3 



which provides the solution 

6a(t) = dtF* + C 2 e- at . 



(14) 



where C\ and Ci are constants of integration, and uj is 
defined by 



po(l +w)(l + 3w) 



(15) 



In order to avoid the blow-up, due to the exponential 
increase in the scale factor, or the collapse, the solution 
is stable in the range 



1< w < -1/3. 



(16) 



Note that this interval violates the strong energy condi- 
tion which stipulates that p + 3p > 0. 

Since the cosmological constant of the classical Ein- 
stein universe is given by 



A 



-p (l + 3to) , 



(17) 



and we are only considering positive energy densities, it 
turns out that it is negative in the region of stability. 

Now, the full modified perturbation differential equa- 
tion, Eq. (O, provides the following solution 



Sa(t) = Cie" l£ + C 2 e- Wlt + C 3 e" 2£ + C 4 e 



(18) 



where Cj (with i — 1...4) are constants, and the param- 
eters oji and u>2 are given by 



Wl 2 



{ 



po(l + w) — 4cra 2 



5 (7 or 



iaa 2 - p (l + w)± 

1/2 



± + «0 EMI + w) + 8aa 2 {3w - 1)]] } , (19) 

respectively. 

In the following, we require the cosmological constant 
to be positive. For a = — 1 and A > no stable solutions 
can be found. Considering a = +1 and A > 0, we have 
the following three stability regions: Firstly, 



AI 



8a 2 < po < 



3(7 + Vl7)a 



8a 2 



Po , ^ 1 

^T7o~ w < 3 





< w < - 
3(23 + V17) 3 



The latter inequality reduces to —0.213 < w < —0.146. 
Finally, 



AIII: po > 



3(7 + V17) 



6<a 2 - Pa , / 1 / 

< w < — —2 

3p 3 \ 



'24a 2 + Po 



po 



(24) 



(25) 



These stability regions are summarized in Table [H and 
depicted in Fig. [1] Note that these results are consis- 
tent with the stability condition fan = cPf/dR 2 > of 
cosmological models at high curvatures [10( • 



B. f(R) oc R+ 1/R theory 



In this section, we use the specific form of 



f(R) = R 



ap 4 1 



R 



-2A. 



(26) 



where a = ±1, and p is considered a positive parameter. 
As in the previous example, we have inserted the factor 
a 2 , which simplifies the calculations and the respective 
notation outlined below. 

The choice of the f(R) given by Eq. (f2"6"| has been 
extensively analyzed in the literature since it was first 
demonstrated to account for the late-time accelerated 
expansion of the Universe without the need for the in- 
troduction of dark energy [3(3]. It has also been used 
in the weak- field limit, as now the 1/R term dominates 
for low curvatures. Unfortunately, it was demonstrated 
that the originally proposed form suffers from instabili- 
ties [31(. However, it was recently shown in [l(| that a 
modification of the sign stabilizes the solution such that 
Irr = d 2 f/dR 2 > 0, as emphasized in the Introduction. 

Considering this case, the unperturbed field equations 
([5]) take the form 



3 ap 

Po = — + -T7T ~ A ' 
ag 12 

1 5CT/J 4 . 
PO = -n ^7— + A . 

a 2 



36 



(27) 
(28) 



which implies that the cosmological constant of the mod- 
ified Einstein static universe is given by 



1 , v CTit 

A= 2^(1 + 3^) + -^- 



(29) 



Applying linear perturbation theory, in the procedure 
outlined in the previous case, we deduce the following 
differential equation 
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Case AI 



Case All 



Case AIII 



8a 2 < po < 



3(7+^17)a z 
2 



Po > 



3(7+yTf)a 2 



., 8 ° ,~ Pn < IW < 4 I -2 
24a^+p — 3 



24c 2 +po 



. 5+3VT7 1 
3(23+yT7) ^ ^ 3 



-2 + 



23+yT7 
7+\/T7 



6a -Po 
3po 



< w < 



TABLE I: Summary of the stability regions in the Einstein static universe for the specific case of f(R) oc R + R 2 theory. 
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FIG. 1: The stability regions, for the case of f(R) oc R + R 2 theory, are depicted between the curves. The thick solid line 
represents general relativity. In the lower left region, which tends to general relativity, the cosmological constant is negative. 
However, in the upper triangular-like shaped region (AI, All & AIII) the Einstein static universe is stable and A is positive, 
compare also with Fig. [2] Note that the equation of state parameter w is strictly negative. 



[2o> 4 - 9(1 + w)(l + 3w) Po ] [o> 4 + 18(1 + w)p ] Sa(t) 

- 36 [2a p 4 - 9(1 + w)p Q ] [op 4 + 18(1 + w)p ] 6a" (t) + 648ap 4 5a<- 4) (t) = . (30) 

I 



In the limit pi —> this differential equations also re- following stability regions 
duces to the general relativistic one, i.e., Eq. (flU|) . Equa- 
tion ([50]) provides the following solution 



6a{t) = Cie" 3 * + C 2 e-" 3t + C 3 e" 4t + C^e"" 4 ' . (31 
The parameters u>z and W4 are given by 
' p 4 a + 18p (l + w) r 4 



^3,4 = 



36 p 4 a 



-{2p 4 a 



-3 



± 



3p (l+w)± 

1/2 



Vpo(l + w) [2 M %(-1 + 3w) + 9p (l + wfi] }} -(32) 



As in the previous example, we verify that for a = — 1 
and a positive cosmological constant, the solutions are 
unstable. Considering a = +1 and A > 0, we have the 



BI 



2p 4 



< Po < 



(5 + v41)m 



2^ 4 - 9po 
3(2/t 4 + 3p ) 

Secondly, 



< w < - 
9 



12 
-6 + 



'3(2^ + 3/70 



po 



(33) 
(34) 



BII : po 

7 + 3\/4T 



(5 + V41> 
12 



1 n 13 + y/U 

< w <-\-2 + \ -j=- 

3(13 + V41) 3 \ V 5 + V4T 



(35) 
(36) 
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Finally 

Bill: p > 

p 4 + 3po 



(5 + V41)m 



12 



9po 



< w <- 




3(2/i 4 + 3 Po ) 



(37) 
(38) 



We summarize these stability regions in Table |TT1 and 
depict the solutions in Fig. [2] Note that, as emphasized 
above, these results are also consistent with the stability 
condition f RR = d 2 f/dR 2 > [13]. 



III. SUMMARY AND DISCUSSION 

The Einstein static universe has recently been revived 
as the asymptotic origin of an emergent universe, namely, 
as an inflationary cosmology without a singularity [24(. 
The role of positive curvature, negligible at late times, 
is crucial in the early universe, as it allows these cos- 
mologies to inflate and later reheat to a hot big-bang 
epoch. An attractive feature of these cosmological mod- 
els is the absence of a singularity, of an 'initial time', of 
the horizon problem, and the quantum regime can even 
be avoided. Furthermore, the Einstein static universe 
was found to be neutrally stable against inhomogeneous 
linear vector and tensor perturbations, and against scalar 
density perturbations provided that the speed of sound 
satisfies c 2 > 1/5 (H)]. Further issues related to the sta- 
bility of the Einstein static universe may be found in Ref. 

M 

In this work we have analyzed linear homogeneous 
scalar perturbations around the Einstein static universe 
in the context of f(R) modified theories of gravity. We 
have considered two specific forms of f(R) and found 
the stability regions of the solutions for the scale fac- 
tor perturbation. The first case considered, namely, 
f(R) (x R+R 2 , was motivated by the fact that, in princi- 



ple, R 2 dominates for high curvatures which is expected 
in the early universe. Secondly, we considered the case of 
f(R) oc R+l/R, which is known to generate a late-time 
accelerated expansion phase [30| . and has been used in 
the weak-field limit, as now the 1/R term dominates for 
low curvatures. The stability regions were parameterized 
by an equation of state parameter w — p/p, and it was 
found that in the context of f(R) modified theories of 
gravity the range of the parameter is greatly enhanced 
relatively to the results obtained in general relativity. 
However, in both cases we analyzed, the equation of state 
parameter was strictly negative, an issue which we hope 
to overcome in future work by considering other mod- 
ified gravity models. These results are consistent with 
the condition fan = d 2 f/dR 2 > 0, for the stability of 
cosmological models [lOj. 

Concluding, we have found that the modified Einstein 
static universe, with a positive cosmological constant and 
matter described by the equation of state, p = wp, can be 
stabilized against homogeneous perturbations, contrary 
to classical general relativity. Therefore, we are lead to 
conclude that, in principle, stable modified gravity so- 
lutions, which are unstable in general relativity [34j . do 
indeed exist. That this is actually possible relies on the 
fact that the perturbations of the metric couple to the 
matter perturbations, see our Eq. (jTTJ) . Similar results 
have been obtained in [29j | where the modified Friedman 
equations in loop quantum cosmology are considered. 
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